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Abstract 

The  performance  of  a  thermoelectric  refrigeration  device,  consisting  of  multi-couple  thermoelec¬ 
tric  elements  and  operating  between  two  heat-reservoirs  at  constant  temperatures,  is  investigated. 
The  influence  of  the  external  and  internal  irreversibilities  of  the  thermoelectric  refrigeration  device 
on  the  performance  of  the  system  is  analyzed.  The  general  expressions  of  the  coefficient  of  perfor¬ 
mance  and  power  input  are  derived  by  introducing  some  dimensionless  parameters  and  variables. 
The  coefficient  of  performance  of  the  refrigeration  device  is  maximized  for  a  given  cooling-load 
and  total  heat-transfer  area  of  the  system,  and  consequently,  the  structure  parameters  of  the  device 
and  the  area  ratio  of  the  heat  exchangers  of  the  system  are  determined  optimally.  The  effects  of  the 
various  parameters  on  the  optimal  performance  of  the  device  are  discussed  further.  The  results 
obtained  here  will  be  useful  for  a  more  detailed  investigation  and  for  the  optimal  design  and  man¬ 
ufacture  of  real  thermoelectric  refrigeration  devices. 

©  2007  Elsevier  Ltd.  All  rights  reserved. 

Keywords:  Thermoelectric  refrigeration  device;  Irreversibility;  Structure  parameters;  Heat-transfer  area;  Perfor¬ 
mance  analysis;  Optimal 


Corresponding  author.  Tel.:  +86  592  2189  426;  fax:  +86  592  2185  119. 
E-mail  address:  jcchen@xmu.edu.cn  (J.  Chen). 

0306-2619/$  -  see  front  matter  ©  2007  Elsevier  Ltd.  All  rights  reserved, 
doi:  10.101 6/j  .apenergy.2007.02.008 


Y.  Pan  et  al.  /  Applied  Energy  84  (2007)  882-892 


883 


1.  Introduction 

Thermoelectric  refrigeration  and  generation  devices  occupy  a  niche  market,  because 
they  are  quiet  and  reliable,  and  friendly  to  our  environment  [1-4].  Nowadays,  the  devices 
have  a  distinct  place  in  medical  applications,  scientific  equipment  and  other  applications, 
where  a  high-precision  temperature-control  is  essential.  However,  the  coefficients  of  per¬ 
formance  of  refrigeration  devices  and  the  efficiency  of  generation  devices  are  low  com¬ 
pared  with  traditional  devices  [5,6].  Therefore,  the  improvement  of  performance  of  the 
thermoelectric  devices  is  an  important  issue  in  the  future  applications  not  only  for  devel¬ 
oping  high  technology  but  also  for  improving  our  quality  of  life. 

Generally,  the  coefficient  of  performance  of  thermoelectric  refrigeration-devices  or  the 
efficiency  of  thermoelectric-generation  devices  is  mainly  determined  by  the  parameter 
Z  =  vl2/(RK)  which  is  called  the  figure  of  merit,  where  a  is  the  Seebeck  coefficient,  R  is 
the  electric  resistance,  and  K  is  thermal  conductance  [7,8].  It  is  well  known  that  the  larger 
the  figure  of  merit,  the  better  the  performance  of  a  thermoelectric  device  and  the  higher  the 
efficiency  of  a  thermoelectric  generator  [1,3,9-12].  In  the  case  of  given  semiconductor 
materials,  the  figure  of  merit  Z  depends  on  the  structure  of  the  device.  It  has  been  proven 

that  when  ln/Sn  =  (/p/ SV)^J Knpp/ (/cppn),  Z  attains  its  maximum  and  the  geometric  config¬ 
uration  of  the  thermoelectric  device  is  optimum  [13],  where  /  and  S  are,  respectively,  the 
length  and  cross-sectional  area  of  the  semiconductor  elements,  p  is  the  electrical  resistivity, 
k  is  the  thermal  conductivity,  and  the  subscripts  n  and  p  designate  the  n-  and  p-type  ele¬ 
ments.  On  the  other  hand,  for  real  thermoelectric  devices,  the  operative  conditions  are  dif¬ 
ferent  from  each  other  and  the  thermal  conductances  between  the  thermoelectric  device 
and  the  external  heat  reservoirs  are  always  finite.  These  questions  are  necessarily  consid¬ 
ered  in  the  optimal  design  of  thermoelectric  devices.  This  has  been  noted  in  recent  years 
[14-17]. 

In  the  present  paper,  we  will  investigate  the  optimal  performance  of  a  thermoelectric 
refrigeration  device  with  external  and  internal  irreversibilities  and  analyze  how  the  differ¬ 
ent  operative  conditions,  the  finite-rate  heat  transfer  between  the  thermoelectric  device  and 
the  external  heat  reservoirs,  the  structure  parameters  of  the  semiconductor  elements  and 
the  area  ratio  of  heat  exchangers  affect  the  performance  of  a  multi-couple  thermoelectric 
refrigeration  device.  The  optimally  operating  regions  of  some  important  performance 
parameters  of  the  thermoelectric  refrigerator  are  determined. 

2.  An  irreversible  thermoelectric  refrigeration  device 

In  thermoelectric  materials,  heat  can  be  transported  or  used  to  generate  electricity 
based  on  the  Peltier  and  Seebeck  effects.  According  to  non-equilibrium  thermodynamics 
[18,19],  when  an  electric  current  with  the  density  J  traverses  a  semiconductor  placed  in 
a  temperature  gradient,  one  has 

V  •  JE  =  -v  •  (icvr)  +  Tj  •  Va  -  J  •  J/a  (1) 

where  JE  is  the  current  density  of  energy  in  the  semiconductor,  T  is  the  temperature  inside 
the  semiconductor,  and  a,  k ,  and  o  are,  respectively,  the  Seebeck  coefficient,  thermal  con¬ 
ductivity,  and  electrical  conductivity  of  the  semiconductor. 

A  multi-couple  thermoelectric  refrigeration  system,  as  shown  in  Fig.  1,  consists  of  a 
semiconductor  thermoelectric  device  and  two  heat  exchangers  at  the  hot  and  cold  sides. 
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Fig.  1.  A  schematic  diagram  of  a  multi-couple  thermoelectric  refrigeration  device. 

The  semiconductor  thermoelectric  device  is  composed  of  many  n-  and  p-type  semiconduc¬ 
tor  legs  that  are  connected  electrically  in  series  by  metal  strips  and  thermally  in  parallel. 
The  thermoelectric  refrigerating-elements  are  assumed  to  be  insulated,  both  electrically 
and  thermally,  from  their  surroundings,  except  at  the  junction-reservoir  contacts.  The  cur¬ 
rent  is  assumed  to  flow  in  one  dimension,  as  along  the  arm  of  the  refrigeration  device.  In 
addition,  it  is  assumed  that  the  Seebeck  coefficients  an  and  ap,  the  electrical  resistivities  pn 
and  pp,  and  the  thermal  conductivities  Kn  and  kp  of  the  n-  and  p-type  semiconductor  mate¬ 
rials  are  independent  of  temperature  in  order  to  obtain  some  analytical  solutions.  When 
the  thermoelectric  refrigeration  device  operates  stably,  using  Eq.  (1)  and  the  boundary 
condition  r1(0)  =  r2(0)  =  ---  =  r/  (())  =  •••  rn(0)  =  rhj  and  r1(/p)  =  r2(/n)  =  ---  =  r, 
(/n)  =  •  •  •  rn(/n)  =  Tc j,  we  can  derive  equations  of  heat  conduction  of  a  multi-couple  ther¬ 
moelectric  refrigeration  device  as  follows: 

Qc  =  cdTq  -  l-I2R  -  K(Thj  -  Tc,)  (2) 

and 

a  =  a/rhj  +  X-i2r  -  ^(rhj  -  rcj)  (3) 

where  Qh  and  Qc  are  the  heat  flows  from  the  device  to  the  heat  sink  and  from  the  cooled 
space  to  the  device,  a  =  (ap  —  otn)n ,  R  =  ( pplp/Sp  +  pJJS^n  and  K  =  (kpSp/Ip  +  KnSjln)n 
are,  respectively,  the  Seebeck  coefficient,  electrical  resistance  and  thermal  conductance  of 
the  semiconductor  couples,  n  is  the  total  number  of  couples,  and  /  is  the  working  electrical 
current. 

According  to  the  important  relation  between  the  structure  parameters  /n/Sn  and  lp/Sp 
mentioned  above,  the  electrical  resistance  and  thermal  conductance  of  a  multi-couple  ther¬ 
moelectric  refrigeration-device  can  be,  respectively,  expressed  as 

R  =  n  (Pp  +  p)  (h/Sp)  (4) 

and 

K  =  n(^Kp  +  '  KpKnpn/ pp^  (iSp/Zp)  (5) 

It  can  be  found  from  Eqs.  (2)  and  (3)  that  the  internal  irreversibility  comes  from  Joule’s 
heat  I2R  due  to  the  electrical  current  and  heat  leak  K(Th LJ-  —  Tc j)  due  to  the  temperature 
difference  between  the  hot  and  the  cold  junctions  of  the  device.  Due  to  the  finite-rate  heat 
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transfer  between  the  thermoelectric  refrigeration-device  and  the  heat  reservoir,  the  temper¬ 
atures  Th j  and  Tc j  of  the  thermoelectric  refrigeration  device  are  different  from  those  of  the 
heat  reservoirs,  Th  and  Tc ,  and  ^  Th  ^  ^  7"cj.  Assuming  that  the  heat  transfer 

between  the  hot  and  cold  junctions  of  the  thermoelectric  refrigerator  and  their  respective 
reservoirs  obeys  Newton’s  law  [18,20],  one  has 

Qc  =  UCAC(TC  -  rcj)  (6) 

and 


Qh  =  UhAh{Thj  -  rh)  (7) 

where  U h  and  Uc  are  the  overall  heat-transfer  coefficients  of  the  heat  exchangers  at  the  hot 
and  cold  sides,  and  Ah  and  Ac  are  the  heat-transfer  areas  of  the  hot  and  cold  junctions. 
By  using  Eqs.  (2),  (3),  (6),  and  (7),  Th j  and  Tc j  can  be,  respectively,  expressed  as 

KjUoAJ^  +  UhAT\)  +  UoAcUhAjTb  +  UhATW  +{K+  \J^jl)Rp_  +  RaP/2 

hJ  ~~  K(UCAC  +  UhAh )  +  UcAcUhAh  -  (a/)2  +  (UhAh  -  UcAc)oc/ 

(8) 

and 

KjUcAcTc  -  UhAiJ\)  +  UbAhUcAjy -  UcAcTcy.I  +  (g  +  Uh,4h/2)f?/2  -  f?a/3/2 

CJ_  K(UCAC  +  U^Ah)  +  UcAcUhAh  —  (a/)2  +  -*  UcAc)cd 

(9) 

Substituting  from  Eqs.  (8)  and  (9)  into  Eqs.  (6)  and  (7),  one  can  obtain 

_  (flh  -  1)(1  -x)x-j(x  +  Ciy)(l  -x)  -|/[l/(zrh)](x  +  Ci7  +  2Ci)(l  -x) 

[(1  —  x)  +  C2  —  C2j]  (x  +  C|/)  —  [(1  —  x)  —  C2j\C] 

(10) 

and 


=  (1  ~  0c)(l  -x)x-7[(l  -x)  -  C2y]x  +  |/[l/(zrh)][(l  -x)  -  C2y  +  2C2]x 

[(1  —  x)  +  C2  —  Cij\(x  +  Ci/)  —  [(1  —  x)  —  C2AC1 

where  qh  =  Qh/(KTC),  qc  =  Qj{KTc),  0h  =  1/0C  =  Th/Tc,  j  =  al/K,  Z  =  y2/(  RK),  C,  =  K/ 
( UhA),  C2  =  K/(UCA),  and  x  =  AjA  are  the  dimensionless  quantities,  while  A  =  Ah  +  Ac 
is  the  total  heat- transfer  area  of  the  heat  exchangers. 

Using  Eqs.  (10)  and  (11),  one  can  derive  the  coefficient  of  performance  and  power  input 
of  the  thermoelectric  refrigeration  system.  They  can  be,  respectively,  expressed  as 

_ (x-  l)x[2Z7h(0h  -j-  1)  +/0b]  +  CJ2[2(x-ZTb)  -  {J-2)(x- l)6h] _ 

C 1  /  (x  -  1 )  [j0b  -  2  ( ZT h  +  9h )]  +  xj 2  8h  [2  +  C2  (2  +/)  -2x]+  IxjZT h  (x  -  1  +  0h  +  C2jdh  -  xdh ) 

(12) 


and 


($h  +  9c)(x—  l)x+/(C,  —  C\x  +  C2x)  +  \j2[\ / (ZT^)\\2x(\  —  x  —  C2  +  Ci)  —xj(C\  +  C2)  +  C  \(J +2)\ 
[( 1  -  x)  +  C2  -  Cj\  (x  +  CJ)  -  [( 1  -  x)  -  Cj\ Ci 


where  p  =  P/(KTC )  =  qh-  qc. 


(13) 


886 


Y.  Pan  et  al.  /  Applied  Energy  84  (2007)  882-892 


3.  Performance  analysis  and  parametric  optimum  design 

For  given  semiconductor  materials  and  specified  operation  conditions,  the  parameters 
ZTh,  UhA,  UCA  and  6h  are  some  known  quantities.  It  is  seen  from  Eqs.  (10)  and  (12)  that, 
if  QjTc  is  given,  the  coefficient  of  performance  is  only  a  function  of  both  K  and  x.  Using 
Eqs.  (10)  and  (12),  one  can  plot  the  three-dimensional  graph  of  the  coefficient  of  perfor¬ 
mance  varying  with  the  parameters  K  and  x,  as  shown  in  Fig.  2,  where  the  parameters 
QJTC  =  0.003  W  K  *,  9h  =  1.2,  ZTh  =  1  and  UCA  =  UhA  =  UA  =  5  W  K  1  are  chosen. 
It  is  seen,  from  Fig.  2,  that  the  coefficient  of  performance  of  the  refrigeration  system  first 
increases  and  then  decreases  as  K  and  x  are  increased.  It  shows  that  there  exist  the  max¬ 
imum  coefficient  of  performance  emax  and  the  corresponding  optimal  values  of  K  and  x  for 
other  given  parameters.  It  can  be  seen  from  Eqs.  (10)  and  (12)  and  Fig.  2  that  for  the  dif¬ 
ferent  choices  of  QjTc ,  K  and  x  will  have  different  optimal  values.  Using  the  extremal  con¬ 
ditions  ds/dK  =  0  and  ds/dx  =  0  and  Eqs.  (10)  and  (12),  one  can  generate  the  smax  ~  QjTc, 
jopt  ~  Qc/Tc,  Kopt~Qc/Tc,  xopt  ~  Qj Tc,  and  emax  -  Kopt  characteristic  curves  of  the 
multi-couple  thermoelectric  refrigeration  system,  as  shown  in  Figs.  3-6,  respectively, 
where  the  parameters  0^  =  1.2  and  ZTh  =  1  are  chosen  and  curves  I  and  II  correspond 
to  the  cases  of  UCA  =  Uh,  A  =  10  W  K-1  and  5  W  K-1,  respectively.  It  is  seen  from  these 
curves  that  the  maximum  coefficient  of  performance  emax  is  a  monotonically  decreasing 
function  of  QjTc ,  while  jQ pt,  Kopt  and  xopt  are  monotonically  increasing  functions  of 
QJTC. 

It  is  seen  from  Eq.  (5)  and  Fig.  4  that,  for  a  given  QjTc ,  the  optimal  values  of  the 
parameters  n(Sp/lp)  and  x  can  be  obtained  from  these  optimal-relation  curves.  For  exam¬ 
ple,  when  6h  =  1.2,  ZTh  =  1,  UCA  =  UhA  =  5  W  K-1  and  QjTc  =  0.003  W  K-1,  the  max¬ 
imum  coefficient  of  performance  occurs  at  K  =  Kopt  =  0.0255  W  K-1  and  x  =  xopt  =  0.62. 
Then,  the  optimum  value  of  n(Sp/lp )  is 

(n^)  = - Q;°255  (14) 

V  w  opt  KP  +  JKpKnpJpp 


Fig.  2.  The  variation  of  the  coefficient  of  performance  £  with  the  thermal  conductance  K  and  area  ratio  x,  where 
the  parameters  QjTc  =  0.003  W  K-1,  0h  =  1.2,  ZTh  =  1  and  UCA  =  UhA  =  UA  =  5  W  K-1  are  chosen. 
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Fig.  3.  The  maximum  coefficient  of  performance  emax  and  the  corresponding  dimensionless  current  jQ pt  versus  Qj 
Tc.  The  values  of  the  parameters  0h  and  ZTh  are  the  same  as  those  used  in  Fig.  2  and  curves  I  and  II  correspond  to 
the  cases  of  UCA  =  UhA  =  10  W  K-1  and  5  W  K-1,  respectively. 


Q/T 

c  c 


Fig.  4.  The  thermal  conductance  Kopt  and  area  ratio  xopt  at  the  maximum  coefficient  of  performance  versus  Qj 
Tc.  The  values  of  the  parameters  0h,  ZTh,  UCA  and  UhA  are  the  same  as  those  used  in  Fig.  3. 

According  to  Eq.  (14),  the  optimal  characteristic  curves  obtained  above,  and  the  techno¬ 
logical  requirements,  one  can  determine  the  optimal  values  of  Sp/lp  and  n.  Consequently, 
Eq.  (14)  and  the  optimal  characteristic  curves  of  the  system  may  provide  theoretical  guid¬ 
ance  for  engineers  to  design  the  structure  of  thermoelectric  devices. 

It  is  also  seen  from  Fig.  4  that  the  optimal  values  of  K ,  Aj A  and  AjA  are  closely 
dependent  on  QjTc.  This  shows  that  for  the  thermoelectric  refrigeration-system  with 
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Fig.  5.  The  maximum  coefficient  of  performance  emax  versus  the  thermal  conductance  Kopt  at  the  maximum 
coefficient  of  performance.  The  values  of  the  parameters  0h,  ZTh,  UCA  and  UhA  are  the  same  as  those  used  in 
Fig.  3. 


Fig.  6.  The  local  maximum  of  the  coefficient  of  performance  £im  versus  the  area  ratio  x.  The  values  of  the 
parameters  6h,  ZTh,  UCA  and  UhA  are  the  same  as  those  used  in  Fig.  2  and  curves  I  and  II  correspond  to  the  cases 
of  QJTC  =  0.003  W  K-1  and  0.004  W  K1,  respectively. 

the  different  operative  conditions  which  are  required  practically,  the  optimal  values  of  the 
structure  parameters  n(Sp/lp)  of  the  thermoelectric  device  and  the  area  ratio  of  the  heat 
exchangers  must  be  different  from  each  other.  This  is  worthwhile  to  consider  when  design¬ 
ing  the  heat  exchangers  of  the  multi-couple  thermoelectric  refrigeration-system. 
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In  addition,  for  given  parameters  x,  QjTc ,  0 h,  ZT h,  UCA  and  U^A,  using  the  extremal 
conditions  (de/dK)x  =  0  and  Eqs.  (10)  and  (12),  one  can  obtain  the  relations  of  the  varia¬ 
tion  of  the  local  maximum  of  the  coefficient  of  performance  and  corresponding  parameters 
K  and  j  with  the  area  ratio  x,  as  shown  in  Figs.  6-8,  respectively,  where  the  parameters 
9h  =  1.2,  UCA  =  U^A  =  5  W  K-1  and  ZTh  =  1  are  chosen  and  curves  I  and  II  correspond 
to  the  cases  of  QjTc  =  0.003  W  K-1  and  0.004  W  K-1,  respectively.  It  is  seen  from  these 
figures  that  the  local  maximum  of  the  coefficient  of  performance  of  the  system  glm  first 
increases  and  then  decreases  while  the  corresponding  optimal  values  j\0  and  KXo  of  the 
dimensionless  current  and  thermal  conductance  first  decrease  and  then  increase  as  x  is 
increased.  It  shows  that  there  is  an  optimal  value  of  x  at  which  8Xm  attains  its  maximum 
value,  while  jXo  and  KXo  attain  their  minimum  values  for  other  given  parameters.  Also 
the  optimal  values  of  x  are  different  for  the  different  parameters  at  the  different  operating 
conditions.  It  can  be  seen  from  the  curves  in  Figs.  6-8  that  there  is  an  important  relation: 

^(Amin)  <  x0pt(smax)  <  x(ymin  )  (15) 

where  xopt(£max)  is  the  area  ratio  at  the  maximum  coefficient  of  performance,  and  x(^min) 
and  x(/m in)  are,  respectively,  the  area  ratios  at  the  minimum  thermal  conductance  and  min¬ 
imum  dimensionless  current.  Figs.  9  and  10  further  show  the  variations  of  8Xm  with  KXo  and 
y'lo,  respectively.  It  is  seen  from  Figs.  9  and  10  that,  when  j\Q  > ymin  or  KXo  >  Km in,  there  are 
two  different  coefficient  of  performance  for  a  given  dimensionless  current  jXo  or  thermal 
conductivity  KXo.  One  always  wants  to  obtain  the  coefficient  of  performance  as  large  as 
possible  for  a  given  cooling  load  QjTc.  When  x  >  x(Xmin),  the  coefficient  of  performance 
of  the  system  is  larger  than  £(Xmin).  Similarly,  when  x  <  x(/min),  the  coefficient  of  perfor¬ 
mance  of  the  system  is  also  larger  than  £(/min).  Consequently,  it  is  found  that  the  range 
of  the  area  ratio  x  should  be 

x(Kmin)  <  x  <  x(Jmin)  (16) 


Fig.  7.  The  dimensionless  current  j\0  at  the  local  maximum  of  the  coefficient  of  performance  versus  the  area  ratio 
x.  The  values  of  the  parameters  0h,  ZTh,  UCA,  UhA  and  QjTc  are  the  same  as  those  used  in  Fig.  6. 
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Fig.  8.  The  thermal  conductance  K\0  at  the  local  maximum  of  the  coefficient  of  performance  versus  the  area  ratio 
x.  The  values  of  the  parameters  0h,  ZTh,  UCA,  UhA  and  QjTc  are  the  same  as  those  used  in  Fig.  6. 


Fig.  9.  The  local  maximum  of  the  coefficient  of  performance  £im  versus  the  thermal  conductance  KXo  at  the  local 
maximum  of  the  coefficient  of  performance.  The  values  of  the  parameters  0h,  ZTh,  UCA,  UhA  and  QjTc  are  the 
same  as  those  used  in  Fig.  6. 


Only  when  x  =  xopt,  can  the  coefficient  of  performance  of  the  system  attain  its  maximum 
and  the  multi-couple  thermoelectric  refrigerator  system  be  operated  in  the  optimal  state. 
In  this  case,  the  thermal  conductance  K,  dimensionless  current  j  and  the  structure  param- 
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Fig.  10.  The  local  maximum  of  the  coefficient  of  performance  eim  versus  the  dimensionless  current  j\Q  at  the  local 
maximum  of  the  coefficient  of  performance.  The  values  of  the  parameters  0h,  ZTh ,  UCA,  UhA  and  QjTc  are  the 
same  as  those  used  in  Fig.  6. 


eters  n(Sp/lp)  of  the  thermoelectric  device  should  be  controlled  to  satisfy  the  following 
conditions: 

K  =  ^lo(-^opt)  =  ^opt 
j  =  7lo(Xopt)  =  iopt 

and 

£p  ^opt  f  Sp\ 

nf  = -  =  (  riy  j  (19) 

P  KP  +  \JKpKnPn/Pp  V  P/<>P‘ 

The  above  results  show  that  gmax,  xopt,  jopt,  Kopt  and  (nSp/lp)opt  are  important  parameters 
of  a  multi-couple  thermoelectric  refrigeration  system,  which  are  very  useful  for  the  optimal 
design  of  real  thermoelectric  refrigeration-systems. 

4.  Conclusions 

The  externally  and  internally  irreversible  multi-couple  thermoelectric  device  model  is 
used  to  analyze  the  performance  of  a  thermoelectric  refrigeration  system.  Some  fundamen¬ 
tal  relations  are  derived,  from  which  several  optimal  performance  characteristic  curves  of 
the  thermoelectric  refrigeration  system  are  plotted.  The  optimal  operating-regions  of  some 
important  performance  parameters  are  specified.  The  optimal  operating-state  of  the  max¬ 
imum  coefficient  of  performance  for  a  thermoelectric  refrigeration  device  has  been  deter¬ 
mined.  The  influence  of  the  thermal  conductance  between  the  device  and  the  heat 


(17) 

(18) 
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reservoirs  and  the  cooling  load  on  the  performance  of  the  system  has  been  analyzed  quan¬ 
titatively.  The  results  obtained  here  reveal  some  general  performance  characteristics  of  real 
multi-couple  thermoelectric  refrigeration  systems  operating  at  various  conditions  and 
could  be  used  to  guide  the  optimal  design  and  manufacture  of  real  thermoelectric- 
refrigerators. 
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